Cross section dependence of event rates at neutrino telescopes 
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We examine the dependence of event rates at neutrino telescopes on the neutrino-nucleon cross 
section for neutrinos with energy above 1 PeV, and contrast the results with those for cosmic ray 
experiments. Simple scaling of the Standard Model cross sections leaves the rate of upward events 
essentially unchanged. Details, such as detector depth and cross section inelasticity, can influence 
rates. Numerical estimates of upward shower, muon and tau event rates in the IceCube detector 
confirm these results. 
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Introduction. Neutrinos with energy above 0(10^) 
GeV have not been detected. They are sought after in 
both neutrino telescopes and cosmic ray experiments for 
two reasons. One, they may point back to (new) sources 
at cosmological distances and reveal the structure of their 
internal engines. Two, they enable the study of neu- 
trino interactions at energies far beyond that achievable 
at colliders. The latter is especially interesting from a 
particle physics point of view because neutrinos only in- 
teract weakly according to the Standard Model (SM). 
Any new physics that affects how particles interact is 
therefore more likely to be detected in the neutrino sec- 
tor. However, practical issues hamper the extraction of 
such information. Knowledge of astrophysical sources 
and fundamental particle interactions must be built up 
simultaneously, since we are largely ignorant of both in 
the energy regime of interest, and knowledge of one is 
dependent on knowledge of the other. 

Currently, the situation is further complicated by con- 
fusion. Cosmic ray experiments detect upcoming leptons 
produced in the Earth by charged current interactions of 
"Earth-skimming" neutrinos. The effective area for de- 
tecting these upward leptons is a narrow projection of 
the actual detector surface area. The cross section de- 
pendence for the rate of such events has been derived in 
Refs. 0,12,13. Subsequently, these results have been em- 
ployed to assess the potential of neutrino telescopes 
where they are inapplicable. Our motivation is to eluci- 
date clearly (via qualitative reasoning, analytical calcula- 
tion and numerical analysis) how event rates at neutrino 
telescopes depend on the cross section and to show why 
the event rates in cosmic ray experiments have a differ- 
ent cross section dependence. We emphasize the role of 
the different cross sections and inelasticity in neutrino 
propagation and in the generation of upward shower and 
upward lepton events in providing probes for new physics. 
We begin by using simple models to determine the gen- 
eral dependence of the event rates on the cross sections 
that govern propagation and detection, and apply these 
concepts to different experimental configurations. We 
then reproduce our main results analytically. Finally, we 
perform numerical calculations of the dependence of rates 
on cross section for the IceCube experiment |6[, which 
confirm the conclusions of the simplified models. 

Qualitative, geometric arguments. Event rates 



are given by r{E^) = ^ = / dflAe^, where is the 
effective area for detection, and we use a diffuse, isotropic 
flux for ^ . We assume that the dimensions of the detec- 
tion volume are small compared to the length scale for 
neutrinos to be absorbed. One may then approximate 
Ae = naiVe, where T4 is the effective volume over which 
an interaction may be detected. Ve can be expressed as 
Ve — Ap£, where Ap is the area of the detector projected 
against the neutrino direction, and £ is the portion of 
the neutrino path to which the detector is sensitive; n 
is the nucleon density within V^', and <7i is the total (t), 
charged current (c) or neutral current (n) cross section 
corresponding to the process which makes the observable 
events. For shower detection, ai = at = Cc + Cn as both 
charged and neutral current interactions produce show- 
ers, whereas for charged leptons ai = Uc- Introducing 
the interaction length Ai = \/{n(Ji), the effective area 
becomes Ae = Ap£/\i = Aplncn. 

We first consider showers from downward neutrinos. 
Since showers are contained, £ is given by the detector 
scale s, so Ae = Aps/Xt- Typically, the detector is at a 
shallow depth d compared to the absorption length for 
neutrinos, so one may integrate dil over the whole sky. 
The event rate for downward showers is 
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which increases with at- The second case is the produc- 
tion of downward leptons. The effective area is given 
by Ae — Ap£/Xc. Here, £ is determined by the lepton 
stopping (dE/dx) or decay, by the chord length to the 
surface, and by the detector size. We define leptons as 
detectable until they have lost 1/e of their initial energy, 
after which they would be difficult to distinguish from 
more numerous lower energy leptons. This is a simple cri- 
terion, which can only be refined with a full experimental 
Monte Carlo. Taking dE/dx = —7 — 13E, the stopping 
range for lepton j is 1/ (3j and 7 characterizes the ioniza- 
tion process. The decay length is Ij = CTjE/rrij, and we 
denote the path length from the detector to the surface 
by l{e). Qualitatively, £ = s/2 + Mm{l/Pj,lj,l{e)), and 
the event rate is approximately 
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which scales as Oc- For simpHcity, we have suppressed 
the dependence on the depth of the detector. 

Next we turn to upward events, and penetration 
through the Earth. We introduce the concept of neu- 
trino attenuation length, Aq, and attenuation cross sec- 
tion a a — ^iVi which is weighted by inelasticity. For 
neutral currents the inelasticity in the energy range under 
consideration, (y) ~ 0.2 — 0.25 0], accounts for down- 
scattering. For Ve and Vy^ charged current interactions 
the charged leptons stop, so y = 1 is appropriate. For Vr 
the interplay between r energy loss and decay gives an 
effective inelasticity in regeneration which increases with 
energy from about 0.6 to 1 For the SM at 10 PeV, 
aa ^ (0.77, 0.77, 0.5)at for {ve, Vp,, Vr). 

It simplifies the discussion to consider the case of a 
shallow detector near the surface for which d = 0. The 
length of the chord traversed by an upward neutrino is 
I = 2Rsm6, where R is the radius of Earth and 9 is the 
entry angle measured from the horizon. The attenuation 
length Aa limits the maximum chord length, so that the 
solid angle over which neutrinos are accepted scales as 
Q = 2TTsm9 = wXa/R- Using Ae — Aps/Xt, the rate for 
upward shower events is 



(3) 



Although the event rate will depend on details of ctq and 
at, to the extent that all cross sections scale in a similar 
manner, the upward shower rate is independent of new 
physics. This is a consequence of the competition be- 
tween two opposing effects: as the cross sections increase 
the number of neutrinos entering the detector from below 
decreases with increased attenuation, but the probability 
for those neutrinos to interact increases. 

The fourth case is upward leptons. We use 
Ae = Api/Xc- Folding in the reduced solid angle. 
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If neutral current and charged current cross sections 
change proportionally, then the upward lepton rate is in- 
dependent of cross section. However, if the new physics 
only adds to the shower rate, then the rate decreases due 
to the decrease in acceptance without a compensatory 
increase in the charged current interaction probability. 

Since it was one of our motivating concerns, we ex- 
amine the rate for cosmic ray air shower experiments to 
detect upward going leptons. In this case the detection 
criterion is that the leptons pass through a horizontal 
planar area A^) ^ corresponding to the footprint of 
the cosmic ray experiment, and the projected area is 
Ap = ylosin6'. Since / = 2i?sin6i < Aa, sin 6* < Xal(2R), 
so the projected area decreases with an additional factor 
of I/ctq. Integrating over solid angle, 
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We thus come to the conclusion of Ref. H that the event 
rate for upward leptons scales as l/cTg, where one factor 
of (Ta comes from geometry, and one factor from ac/<Jg 
assuming that the new physics does not increase the 
charged lepton production. However, we also see that 
this result does not apply to an in-ice detector such as 
IceCube, where the event rate for shower and charged 
lepton detection scales as at/cTg and ac/^g respectively. 

An oversimplification that captures the essential point 
of the above analysis for upward leptons is as follows. 
Suppose a detector can be described as a box with height 
h and side s. The horizontal surface area is Aq = and 
the effective geometric area for incident direction 9 is 
Ap ~ sh cos 9 -h Aosin6' ~ s{h + s9). A detector like 
IceCube has h s and Ap ~ which is independent 
of 9. Cosmic ray detectors have a planar geometry for 
which h <^ s9 and hence Ap ~ Aq9. It is this additional 
factor of 9 that yields the extra factor of l/tTa in Eq. ||SJ). 

Analytic model. The probability that an upcom- 
ing neutrino propagates a distance x along a chord 



is P{x) 



and the probability that it inter- 



acts and produces a detectable signature i in the 
interval dx is dxjXi. The rate of interactions is 
Y^Ap% //^ 27r^ e"^/^»f2^ ^hich on integration 
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limit that 2i? > Aa > ^, F ~ t^^p^-Ii^' which agrees 
with Eqs. © and (gl). 

For a surface detector an additional factor of sin 9 — 
1/{2R) must be included in the dl integral to project out 
the normal component of the lepton flux emerging from 
the Earth. For 2i? > Aa > we recover Eq. ((S)). 

Numerical results. So far, we have discussed neu- 
trino event rates in various settings, but with several 
simplifying assumptions. We have used a shallow de- 
tector, a constant density Earth model, and have sim- 
plified the treatment of regenerated and downscattered 
neutrinos. To test that our conclusions are robust, we 
perform numerical integrations to estimate the event 
rates for upward showers and leptons produced by an 
isotropic total flux of neutrinos and antineutrinos ^ = 
6 • 10-8(i;^/GeV)-2 (cm^.s.sr.GeV)-! || with an as- 
sumed flavor ratio at Earth of 1:1:1. The initial neutrino 
energies extend from 10^ to 10^^ GeV. This model, used 
in Refs. for illustration only. A well- motivated 

alternative model produces a considerably larger flux rel- 
evant to IceCube [lO| . 

We numerically solve the coupled Boltzmann equa- 
tions TT] for three neutrino species and r leptons, using 
a realistic Earth model to determine the flux of neutri- 
nos at the detector. We include downscattering, dE/dx 
for r's, and neutrino production in r decay. For numer- 
ical integration we divide the flux into 120 logarithmic 
energy bins. We solve for the flux in each of 41 angular 
bins, logarithmically spanning a range of 9 from ~ 0.1 
to 90 degree. The angular binning allows us to focus on 
narrowing angular ranges near the horizon as the cross 
sections increase. 
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Event rates are determined using ice as the detector 
medium for both showers and propagating charged lep- 
tons. We do not attempt to distinguish different event 
topologies, e.g., double bang, lollipop. Rather, we keep 
the simple criteria of the previous sections, i.e., that all 
showers are detected and that leptons have a practical 
range ^ = s/2 + Min(l//3j, /j). We use Ap = 1.4 km^ 
and Ve — 2.25 km'^, IceCube values appropriate to 
the middle of the energy range under consideration, and 
infer s = 1.6 km. We use = 3.9 • 10^^ cm^^ and 
Pr = 0.78 • 10~^ cm~^, average values for ice in the same 
energy range. 

We integrate across angles and present rates in log- 
arithmic energy bins from E^, = 10^-^ — 10^ GeV. The 
Alogif^ El, — 0.5 bin width is motivated by the even- 
tual need to separate a smoothly changing flux and SM 
cross section from a hypothetical new physics cross sec- 
tion which may 'turn on' over approximately a decade in 
energy. Our energy bins are for neutrino energy at the 
detector. A given bin may include down scattered neu- 
trinos, or neutrinos produced in propagation by r decay. 

We parameterize physics models by three numbers ac, 
an , and . The first two reflect the strength of charged 
and neutral current reactions, at = ataf^, normalized to 
the SM total cross section. For our calculations, we take 
the inelasticity dai/dy to be the same as in the SM, and 
assume that the Ui are independent of energy. For a'^ 
we deflne a new neutral current process also normalized 
to crf'^, but which is completely inelastic. Thus, the SM 
would be described by (ac,Q;„,a^) = (rc,r„,0), where 

In Fig. n] we show the expected number of IceCube 
events per year as a function of a for case I: a model 
with (ac,a„,Q;J^) — (arc, ar„, 0), such as might occur 
if QCD saturation effects alter t he g rowth predicted for 
neutrino nucleon cross sections tl2|. The upper panel 
shows shower rates in the five energy bins as a function 
of cross section for the nominal IceCube depth of 1.9 km. 
For moderate values of a the rates are essentially inde- 
pendent of cross section, confirming Eq. (jS)). For low 
values of a the rates decrease because the Earth is no 
longer opaque, and at some point there is no compen- 
sating increase in fiux to account for the reduced cross 
section. For large values of a the rates de creas e because 
the horizontal distance to the surface, V^dR, becomes 
comparable to the attenuation length Xa. To confirm 
this, we redid the calculations, except now with a shal- 
low detector {d — 0) just below the surface. The results 
arc consistent with a fiat event rate, even for large a. The 
slight increase in rates at large a is due to a decreased 
Earth density (and increase in Xa) sampled by the neu- 
trinos as their paths are constrained to lie closer to the 
surface. Since the rates are independent of cross section, 
the integrated rate in a bin drops as l/E^, for an i?^^ 
differential spectrum. 

The lower panel shows similar results for fi and r lep- 
tons, confirming Eq. The main feature to note is that 
at low energy the r rates are low since r decay limits £ to 
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FIG. 1: Expected event rates at IceCube for lO'''^ < < 
10^ GeV (in five bins of width A log^Q Ei, = 0.5) for variation 
in strength of SM cross sections. Top: showers, solid - 1.9 km 
deep detector, dot - shallow detector. Bottom: leptons, solid 
- /i, dot - r. In each case, the uppermost curve corresponds to 
the lowest energy bin. We have adopted the flux of Ref. 



the detector size, while at high energies r dominates since 
decays are delayed and (3t < Different assumptions 
about event identification {e.g., requiring double bang for 
t) would lead to different relative rates for /j, and r. Inte- 
grating our results above 1 PeV, we find an expectation 
of 2.4, 2.2 and 0.9 events per year for upward shower, 
muons and taus, respectively. 

In Fig. 121 we show results for case II, a scenario where 
the SM cross sections are unaffected, but new interac- 
tions produce showers, i.e., {ac,an,an) = ifr., Tr,., a). An 
example would be a low scale gravity model |l3j | where 
microscopic black holes are produced by neutrino- nucleon 
interactions and then decay 

As for case I, Fig. 12 shows a shower rate that is es- 
sentially flat for moderate values of a. The rate for lep- 
tons, however, falls inversely with a because ac is fixed 
at the SM value while CTq is changed by new physics; see 
Eq. Q . The neutrino fiuxes decrease because absorption 
increases, but this is not compensated by an increasing 
charged current interaction rate. 

Figure [3 also show results for a model with 
(ttc, a„, aj^) = {rc,rn + a,0), i.e., where new neutral cur- 
rent interactions are added with inelasticity as in the SM, 
such as might occur due to graviton exchange. The dif- 
ference between the two scenarios is due to differences 
in downscattering. In the first case neutrinos are com- 
pletely absorbed during propagation, while in the second 
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FIG. 2: Expected event rates for lO''-* < < 10^ GeV at 
IceCube for the case where new neutral current processes are 
added. For the "black hole" case (lower of each pair) the 
shower rate (solid) is independent of cross section while the 
rates for leptons (dash - /i, dot - r) fall with attenuation. For 
the case of "graviton exchange" (upper of each pair) reduced 
inelasticity and attenuation lead to higher rates. 

they typically lose just 24% of their energy (at 10'' GeV) 
per interaction. Define the normalized attenuation cross 
section as Wa = (Ta/cf^. Then the two scenarios give 
w'g^ ~ 0.77 + a, and Wa — 0.77 + 0.24a. For example, 
consider a = 9 on the right hand edge of Fig. |21 In this 
case w'^/wa — 3, a ratio which is confirmed by the pairs 
of model curves in the figure. 

Summary and discussion. We have made sim- 
ple scaling arguments to describe the expected detection 
rates for high energy neutrinos, with a focus on the Ice- 
Cube experiment geometry. To validate our approach, we 
have performed numerical calculations which propagate 
a diffuse neutrino flux through the Earth. 

Our first result is that the event rates for upgoing neu- 
trinos in volume detectors are independent of the overall 
strength of the cross section. This conclusion is valid for 
over a decade change in the cross section normalization. 

Our second conclusion is that the details matter. The 
event rates depend on a competition between neutrino 
attenuation and interaction lengths. The attenuation 
length depends on inelasticity y as well as cross section, 
so models which alter y from the SM value will also al- 
ter the scaling conclusion. In the extreme case of new 



totally inelastic showering interactions, the lepton rates 
are suppressed by a factor ^ 1/(1 -I- a) as can be seen 
from Fig. H 

Our third conclusion is that the definition of an "event" 
can lead to apples versus oranges comparisons. In our 
analysis, we required muons to retain most of their pro- 
duction energy, but did not require tau leptons to have 
an identifiable tag. On inclusion of tau event selection 
criteria and the implementation of tau energy losses as 
in Ref . fT3| , the expected SM rate for upward going taus 
with energy above 1 PeV becomes 0.3 events/year. This 
is significantly smaller than the expectation in Refs. ^|^, 
whose extrapolation of the results of Ref. 0| above 
10* GeV to values down to 10^ GeV overestimates the 
tau flux at the detector by an order of magnitude. 

The focus of this research is to understand the possi- 
bilities for inferring neutrino-nucleon cross sections when 
the neutrino flux is itself unknown. By taking a ratio of 
two separate measurements, a diffuse flux can be elimi- 
nated. If the ratio differs from the SM prediction, then 
something new is going on. The converse is not true 
- varying the SM cross sections, for example, would not 
affect the ratio of upward showers to upward muons. Sev- 
eral remedies suggest themselves, a) Measure downward 
leptons; however, bundles of muons from cosmic ray air 
showers may provide an annoying background, b) Mea- 
sure the lepton rate near the horizon. The angular distri- 
bution is squeezed to the horizon as a increases, c) Con- 
sider downward showers, which by our estimates form the 
majority of events. Perhaps the cleanest result would be 
a measurment of upward leptons with good angular reso- 
lution near the horizon, normalized to a total number of 
showers dominated by downward events. IceCube seems 
well designed to accomplish this, but the feasiblity of such 
methods is an experimental question, made difficult by 
an unknown flux and anticipated low statistics. 
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